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Abstract 

The anisotropy parameter of two-dimensional equilibrium clusters of site perco- 
lation process in long-range self-affine correlated structures are studied numerically. 
We use a fractional Brownian Motion(FBM) statistic to produce both persistent and 
anti-persistent long-range correlations in 2-D models. It is seen that self affinity makes 
the shape of percolation clusters slightly more isotropic. Moreover, we find that the 
sign of correction to scaling term is determined by the nature of correlation. For per- 
sistent correlation the correction to scaling term adds a negative contribution to the 
anisotropy of percolation clusters, while for the anti-persistent case it is positive. 

1 introduction 

The shape of clusters in percolation type models is an interesting physical property. Numer- 
ical simulations on random models has revealed that percolation clusters at the percolation 
threshold are asymmetric [H El El IH [5] which is a significant result, especially because of 
the close connection between percolation model and the theory of second-order phase tran- 
sition [HE]. 

Generally speaking, the shape of a 2-dimensional cluster is determined by -R^j„ and 
-^maa; i^rnin — ^max) ^ eigenvalues (the principal radii of gyration) of the cluster radius 
of gyration tensor. If -R^j„ = Rmax^ the cluster is spherically symmetric. Otherwise, it is 
anisotropic and we can probe the degree of its anisotropy by defining a proper anisotropy 
quantifier. Family et al pl|, proposed an convenient asymmetry measure, Ag = Rmin/ ^maxy 
called the anisotropy parameter of a s-site cluster. The quantity Ag when properly averaged 
over all clusters with the same size is denoted by (Ag) and is an estimate of the anisotropy 
parameter of s-site clusters in the ensemble. The case (Ag) = 1, corresponds to spherical 
symmetry. For anisotropic objects, (Ag) is less than unity (the term anisotropy parameter 
may be misleading; the shape of the cluster is more isotropic for larger value of Ag). The 
asymptotic behaviour of (Aao) is obtained by taking the limit s oo. Using this method for 
two dimensional random percolation, they observed for the first time that percolation clusters 
are not isotropic and estimated (Aao) = 0.4 as the asymptotic value for the anisotropy of 
infinitely large percolation clusters. 

The most studied percolation systems have been involved with totally random and un- 
correlated models. But, in some important practical applications of percolation theory, the 
nature of disorder is not completely random and there are correlations in the properties of the 
medium. For example, percolation theory has been used to understanding and explaining 
some important aspects of multi-phase flow in porous media [7]. However, studies have 
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demonstrated the existence of long-rang correlations in the permeability distributions and 
porosity logs of some field scale natural porous media like sedimentary rocks [H [9l [TO] . 

A common way of introducing long-range correlations into a medium property is fractional 
Brownian motion (FBM) Db{'^) pTl IT2]. FBM is an ergodic, non-stationary stochastic 
process whose increments are statistically self-similar such that its mean square fluctuation 
is proportional to an arbitrary power of the spatial displacement it 

(DB(l^)-DB(0)]2)~|irp^ (1) 

H is called the Hurst exponent and determines the type of correlations. If if = 0.5, the 
above equation produces the ordinary Brownian motion, which means that in this case there 
is no correlation between different increments. If if > 0.5, then FBM generates positive 
correlations, i.e. all the points in a neighborhood of a given point obey more or less the same 
trend. If ii < 0.5, FBM is anti-persistence, i.e. a trend at a point will not be likely followed 
in its immediate neighborhood. 

In this paper we address the shape of 2-D percolation clusters of self-affine models by 
evaluating their anisotropy parameters. The reason that we have chosen FBM process is 
tenfold. First, FBM generates long-range and at the same time isotropic correlations in 
the field. Therefore, the host lattice retains its isotropy. Second it has been demonstrated 
that such process has practical applications in earth sciences and also reservoir engineering, 
where the permeability field and also the porosity distribution of many real oil reservoirs 
and aquifer follow FBM statistic [131 HI]- A FIB has been used by Schmittbuhl et al [ISJ 
and Sahimi pjlj for generating a percolation model with long-range correlations. 

The paper is organized as follow. After describing the simulation method in section 
2, we present and discuss the results of our extensive numerical simulations of percolation 
processes for both random and self affine models in section 3, with concluding remarks at 
section 4. 



2 Simulation method 

We start with a. L x L square lattice, and assign to each lattice site xt a random number 
drawn from a normalized FBM distribution. There are a number of methods which are 
capable of producing the FBM statistics [131 [12]. We have used one of the most popular 
one, the method of fast Fourier transformation (FFT) filtering which is based on the fact 
that the power spectrum of FBM is given by: 

where oq is a numerical constant, 00 = {001,002), with ooi being the Forger component in the 
ith direction and p = H + 1. In the FAT method, one starts with a white noise W{x,y) 
defined on the lattice sites. The power spectrum of W{x,y) is constant and independent 
of frequency. Therefore, filtering W{x,y) with a transfer function JS{u) generates another 
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noise whose spectral density is proportional to S{uj). The method is straightforward and 
fast, but it usually produces periodic noises. Therefore, one has to produce a larger lattice 
and keep only a portion (typically 1/4 in two-dimensional lattices). In fig.l two different 
realizations of FIB on a 256 x 256 square lattice have been shown. 

To generate percolation clusters one chooses a threshold r. All lattice sites whose 
values are less than r are considered filled. The concentration p of filled sites depends on 
both 0^(1?) and r through: 

rr 

p = T>Bi'^)dlif (3) 
Jo 

The percolation threshold Pc has been defined as the minimum concentration p of filled 
sites to form an infinite cluster with probability 1. It has been observed that percolation 
transitions on long-range self-affine structures behaves differently from ordinary percolation 
in some important aspects [15]. For example, Marrink et al [16] found that percolation 
thresholds Pc of self-affine lattices are strongly dependent on the spanning rule employed 
even when the lattice is effectively infinite. Following the work of Reynolds et al [IT] on 
un-correlated square lattices, they considered three different spanning rules for percolation 
on self-affine models, TZq , the probability of spanning either horizontally or vertically or 
both, TZi, the probability of spanning in a specified direction (e.g. horizontally), and 7^2, 
the probability of spanning both horizontally and vertically. The numerical results obtained 
by Marrink et al established that unlike the urn-correlated random model three critical 
concentrations: p^**^, p^P'^'^, and pl°^^ (corresponding to percolation rules TZq, TZi, and 7^2 
respectively) do not converge for infinite self-affine models. Therefore, the definition of 
the percolation threshold depends on the desired applications and one must consider the 
appropriate percolation threshold to obtain critical parameters of percolation transition of 
self afl&ne models. 

At each specified percolation threshold, there exists a distribution of distinct, finite clus- 
ters (equilibrium percolation clusters). In this work we made the search for clusters by using 
the Hoshen-Kopelman algorithm [TS]. To remove the effect of boundary conditions from our 
numerical estimations, we only count those clusters which do not touch lattice boundaries. 
We call them internal clusters. For each internal cluster of an arbitrary size s, we evaluate 
the cluster radius of gyration tensor G 

N 

G = ^(a^,'l-ir,ir,) (4) 

i=l 

In the above definition, iti is the distance of occupied site i from the cluster's center of mass 
and s is the size of the cluster. The principal radii of gyration of the cluster, -R^j„ and R^axi 
are obtained via diagonalization of G. At the percolation threshold the variations in the i?^, 
have the following asymptotic form: 

(i?^) = riS^^{l + aiS~^ + hiS~^ + ...) i = rain or max (5) 

where v is the leading scaling exponent and 6 is the non-analytical correction-to-scaling 
exponent of the cluster. The coefficients rj, Oj, and hi are all independent of s [1|. The shape 
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of the cluster is then characterized by evaluating its anisotropy parameter, as described 
previously. Finally, we obtain the mean values of anisotropy quantifier by averaging it over 
all the percolation clusters with the same size s, resulted in different lattice realizations. 
Then, the results have been lumped together at the centers of blocks of size [2™~^,2"^]. 
This procedure not only helps to eliminate correction-to-scaling for small clusters [1], but it 
produces new data points which are usually less correlated than the original data [T9] . 

The most time-consuming step in these kind of simulations is generation of FBM statis- 
tics itself. Generating FBM on large lattices requires significant CPU time and computer 
memory. On the other hand, in order to eliminate finite-size effects from numerical results, 
we need lattices with large L . In this work, we have fixed the linear size of our square 
lattices to L = 256 which is both computationally tractable and at the same time large 
enough to make sure that our results are not affected by finite-size effects [TO]. To achieve 
highly accurate results, we estimate the mean values by sampling a large number of lattice 
realizations (30 000 samples in each case). 

In order to make an appropriate and reliable sampling of cluster shapes in configurations 
space it is necessary that the medium linear size be large enough such that all possible 
configurations, including the most anisotropic ones potentially can happen (the condition of 
effectively infinite large medium). This puts a maximum on cluster size s. Indeed, for large 
clusters , i.e. when s >> L"^, the condition of effectively infinite large medium would not be 
fulfilled and as such, the sampling would be in favor of more isotropic configurations. 

3 Results 

Let us first look at the shape of percolation clusters in completely random media. As 
mentioned in the first section, the estimated value of A^o at percolation threshold has been 
reported by several authors. Here, we provide our estimations of the anisotropy parameters 
of percolation clusters for some selected values of p in Tablel. The total number of lattice 
realizations has been 100 000 in this case. The cluster numbers is a function of p, but it 
always decreases rapidly with s [20j. Therefore, when |p — Pc| is not small, the total number 
of clusters of larger sizes is not large enough to give a reliable statistics and we do not include 
the anisotropy parameters of such clusters in our reports. 

As seen from this results, in ordinary percolation the cluster mean anisotropy parameters 
is a function of both p and s. This is natural as the number of a specific cluster configuration 
in random percolation model is proportional to p*(l — p)*, where s is the cluster size and t 
is its perimeter [20|. In fact more careful investigations show some important trends in the 
behaviour of {As{p)) as a function of s and p. For example, looking at the table we can 
see that all the values of s, the mean anisotropy parameter grows with p. Meanwhile, the 
variation of (Ag) for a given p is much more complicated. It seems that for p < p^, (Ag) 
has a maximum at some moderate value of s, while when p > pc, the larger clusters are 
always more symmetric than smaller ones. The behaviour of (Ag) at percolation threshold 
p ~ 0.593 is consistent with equation2, i.e., it evolves to an asymptotic value at large values 
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Table l.The anisotropy parameters of random model 

~m p = 0.454 p = 0.550 p = 0.593 p = 0.650 
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.340(4) 


.355(4) 


.362(5) 


.375(6) 
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.3307(2) 


.3467(2) 


.3545(3) 


.3685(5) 
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.3319(2) 


.3511(4) 


.3607(6) 


.3786(9) 
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.3330(1) 


.3565(3) 


.3688(4) 


.3914(7) 
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.3313(2) 


.3594(1) 


.3752(3) 


.405(1) 
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.326(1) 


.3591(1) 


.3794(2) 


.419(3) 
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.3569(2) 


.3813(3) 





of s [DElEliaiS]. 

To produce the percolation clusters of the self-afiine model we have used the percolation 
thresholds provided in Ref. [16]. In Fig.2 we have shown our estimations of the concentrations 
p of filled sites as a function of cutoff r for three different values of Hurts exponent, H = 0.2, 
H = 0.5 and H = 0.8. For random model, p = r, but as it is seen this is not the case for 
self-affine models. In fact, the fraction of occupied sites at both pl^*^^, and p^J"^'^ in self-affine 
models is less than the corresponding values of random model. This is because compared to 
the random case, the existence of self-affine long-range correlations is in favor of formation 
of larger clusters. Calculation of A^^ = J2k>s''^k, where Uk is the total number of clusters 
with size k per lattice sites (fig. 3) for the case of percolation on either direction shows this 
phenomenon clearly. For comparison, the A^^'s in random model when the concentration of 
filled sites is equal to p^**'* of self affine models with H = 0.2, H = 0.5 and H = 0.8 are 
shown too. 

Our numerical results for the values of the mean anisotropy parameters at percolation 
transition in either direction (rule TZq) and percolation transition in specified direction (rule 
TZi) have been appeared in fig. 4 and fig. 5 respectively. For comparison the corresponding 
values for random percolation have been also included. Interestingly, the differences between 
these two groups of curves are almost negligible, although the difference between pf^'^ and 
pspec jg significant. In fact more investigations suggest that unlike random model, the shape 
of clusters of the self-affine model do not vary significantly with concentration p of the 
filled sites. We may conclude that the introduction of isotropic self-affinity subsidizes (and 
even maybe removes) the dependence of the number of a specific cluster configuration to 
concentration of filled sites. 

One can observe that another important effect of self-affinity is to increase the isotropy 
of percolation clusters. Meanwhile, the difference between the anisotropy parameters of 
different if's and random percolation clusters are not significant especially when H ^ 0. 
This result is in agreement with the previous findings which say that in two dimensions, the 
scaling properties of self-affine model does not differ so much from random model and the 
difference between the critical exponents of two models vanishes as if — > [TU] . 

While the presence of both kind of long-range correlations increases the isotropy of self 
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Figure 2: The concentration p of filled sites as a function of cutoff r for random model (black 
curve), and three different self-affine models: H = 0.2 (green), H = 0.5 (blue) and H = 0.8 (red). 
In each case, the value of r^**^ at p^**^ has been marked (circles). 
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Figure 3: Variation of Ns = J2k>s''^k, where is the normalized cluster number of size k for 
percolation in random model (black) and percolation on either direction in three different self- 
affine models: H = 0.2 (green), H = 0.5 (blue) and H = 0.8 (red). The dash-dot curves are the 
Ng^s in random model when the concentration of filled sites is equal to the pf'^^ of self affine model 
when H = 0.2 (green), H = 0.5 (blue) and H = 0.8 (red). In all cases we have only counted the 
internal percolation clusters of 256 x 256 square lattices. 
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Figure 4: The estimated values of mean anisotropy parameters at p^**'^ for self-affine models with 
H = 0.2 (green), H = 0.5 (blue) and H = 0.8 (red). The results have been lumped together at the 
block centers of size [2™"^, 2™]. For comparison the corresponding values for random percolation 
at the percolation threshold have been also included (black). 
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Figure 5: The estimated values of mean anisotropy parameters at p^^^^ for self affine models with 
H = 0.2 (green), H = 0.5 (blue) and H = 0.8 (red). The results have been lumped together at the 
block centers of size [2™'~^,2™]. For comparison the corresponding values for random percolation 
at the percolation threshold have been also included (black). 
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affine percolation clusters, there is still a significant qualitative difference between the effect of 
persistent and anti-persistent correlations on the shape of clusters. Considering the relation: 

i-^s) — ^oo[l + {0"inin ~ O'max)^ ^ + (&mm ~ bmax)^ ^] (6) 

we can deduce from the behaviour of anisotropy parameters that the sign of correction to 
scaling term, {amin — 0'max)s~^, depends on H: it is positive for H < 0.5 and negative 
for H > 0.5. At H = 0.5 this term is almost zero. Hence, for anti-persisting long-range 
correlations the correction to scaling coefficient is larger than a^ax, which means in 

this case grows faster than R^^x- Again, we see that anti-persistent self-affine models 

shows more similarity with random model (see fig. 4 or fig. 5). In the presence of persisting 
correlations R^ax grows faster than and the anisotropy parameter decreases with s. 

Naturally, for the case H — which separates these two kinds of correlations, the growth 
rate of and R^ax should be equal. 

4 Concluding remarks 

We calculated numerically the anisotropy in the shape of equilibrium percolation clusters in 
long range self-affine correlated square lattices and observed some interesting results. We 
saw that compared to random percolation, the shape of percolation clusters become slightly 
more symmetric and the exact value of anisotropy parameter is a function of Hurts exponent. 
On the other hand, the type of correlations in the medium property, determines the sign of 
non-analytical correction to scaling contribution into the anisotropy parameters of clusters. 
In accordance with the previous studies we also observed that the percolation clusters of 
FBM models become much more similar to those of random models as — > 0. 

Finally, we mention again that, the size of largest cluster that we may consider in our 
analysis is restricted by the lattice size. To estimate the accurate values of or non- 
analytical correction to scaling exponents one should extend the simulation method to larger 
values of L. 

The author would like to thank M. Sahimi for useful hints. 
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